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摘要 : 为 了 获得 大 视 场 的 高 精度 天 文 图 像 ， 实 施 了 一 种 望远镜 CCD 图 像 的 拼接 方法 。 从 
原始 图 像 到 最 终 合成 图 像 的 坐标 转换 采用 了 结合 星 表 UCAC4 的 六 常数 模型 。 不 同 于 硬件 级 的 
拼接 ， 使 用 逐个 像素 的 灰 度 值 再 分 配 的 方案 进行 图 像 融合 。 进 一 步 采 用 云南 天 文 台 24m BR 
镜 拍摄 的 CCD 图 像 进行 了 试验 。 结 果 表 明 ， 该 算法 可 以 产生 较 高 质量 的 大 视 场 CCD 图 像 ， 可 
以 直观 地 发 现 运 动 目标 ， 上 暗 星 信 嗓 比 有 显著 改善 。 高 精度 的 图 像 拼 接 还 与 原始 数据 扭曲 改正 的 
预 处 理 密 不 可 分 。 与 未 做 扭曲 改正 相 比 ， 图 像 拼接 的 位 置 精度 提高 了 约 一 倍 ( 约 0. 02 pixel) 。 
关键 词 : 拼接 ; 大 视 场 ; 几何 扭曲 ; CCD 图 像 
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长 焦距 望远镜 观测 得 到 的 小 视 场 (如 5' x 5) CCD 图 像 常 常 出 现 目标 区 域 参 考 星 不 足 的 问题 '"， 
这 是 由 参考 星 太 暗 或 参考 星 数 太 少 造 成 的 。 文 [2-3] 用 重 迭 露 光 观测 获得 一 系列 CCD 图 像 ， 然 后 将 
这 些 来 自 不 同 平面 的 图 像 归 算 到 同一 天 球面 ， 解 决 了 小 视 场 天 体 测 量 定 标的 问题 。 类 似 地 ， 局 部 平 差 
方法 也 能 较 好 地 解决 这 一 问题 ”。 但 是 这 些 方法 的 归 算 过 程 比较 复杂 。 本 文 直 接 用 CCD 图 像 拼 接 得 
到 大 视 场 的 图 像 并 能 用 于 高 精度 的 位 置 测量 。 

借鉴 一 种 称 作 Drizzle 的 图 像 拼 接 方 法 。 这 种 方法 最 先 由 文 [5 ] 提 出 并 用 于 拼接 Hubble Deep Field- 
North( HDF-N) 的 抖动 观测 图 片 ， 也 用 于 其 他 望远镜 拍摄 的 CCD 图 像 '“ 。 Drizzle 使 用 平移 相 加 能 确保 
高 分 辩 率 ”。 也 有 其 他 一 些 线性 拼接 图 像 的 方法 ， 如 文 [8] 的 方法 ， 文 [9] 的 Imcom 方法 , 文 [10] 的 
iDrizzle 方法 以 及 文 [11 | 的 直接 模型 拟 合 拼接 方法 。 

哈 支 空间 望远镜 ( Hubble Space Telescope, HST) 和 Spitzer 空间 望远镜 (Spitzer Space Telescope, 
SST) 的 图 像 拼 接 可 用 于 解决 原始 图 像 的 欠 采 样 问 题 ”。 地 基 望 远 镜 受 大 气 宁静 度 的 影响 ， 基 本 不 会 
欠 有 来 样 。 不 需要 通过 亚 像 双 级 的 拌 动 观测 ， 只 需 采 用 简单 的 平移 和 灰 度 值 相 加 ， 就 能 保持 拼接 图 像 的 
高 保 真 和 光学 统计 '""。 如 文 [8] 所 述 ， 由 于 光学 设备 的 设计 ( 如 透镜 受 重力 的 影响 ) 会 造成 图 像 的 扭 
曲 ， 所 以 各 个 位 置 的 灰 度 值 未 必 是 其 真实 灰 度 值 。 因 而 对 图 像 做 扭曲 改正 是 必要 的 。 本 文采 用 文 
[12] 的 几何 扭曲 求解 方法 得 到 无 扭曲 图 像 。 该 方法 已 经 成 功 应 用 于 行星 卫星 (Phoebe) 和 近 地 小 行星 
( Apophis ) 的 高 精度 位 置 测 量 ， 效果 显著 :2 4 。 

拼接 图 像 有 如 下 优点 。 第 一 ， 拼 接 图 像 能 用 于 探测 暗 弱 星 像 。 在 拼接 的 图 像 中 ， 同 一 位 置 随 着 重 
合 次 数 的 增加 (相当 于 累积 曝光 时 间 的 增加 )， 信 品 比 不 断 提高 ; 第 二 ， 多 帧 短 曝光 图 像 拼 接 可 以 克 
服 长 时 间 曝 光 的 局 限 ， 确 保 视 场 亮 星 不 饱和 溢出 一 ; 第 三 ， 能 直观 地 发 现 运 动 目标 ; 最 后 还 可 以 
尝试 将 拼接 图 像 用 于 高 精度 位 置 测量 ， 弥 补 小 视 场 参考 星 太 少 的 缺陷 。 

本 文 第 1 市 是 算法 的 描述 ; 第 2 市 对 观测 资料 和 预 处 理 进行 说 明 ; 第 3 世 对 拼接 图 像 的 输出 结 
进行 分 析 ; 最 后 一 节 是 结论 。 
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1 A 法 


在 Windows/ Visual C++ 环 境 下 设计 了 一 球 软 件 实现 相关 算法 。 输 入 是 多 帆 FITS 图 像 及 每 幅 图 中 
主要 参考 星 ( 例 如 UCAC4 星 表 … ) 信息 。 星 像 与 星 表 星 的 匹配 采用 了 文 [18] 的 方法 ， 对 匹配 结果 进 
行 归 算 ， 求 出 扭曲 和 矩阵。 根据 扭曲 矩阵 对 输入 图 像 做 扭曲 改正 得 到 不 扭曲 的 图 像 ， 方 法 详 见 文 [12， 
19] 。 如 果 有 足够 的 参考 星 ， 可 以 使 用 高 阶 常数 模型 进行 位 置 归 算 。 但 实际 上 ， 只 要 扭曲 模型 稳定 ， 
可 以 由 密集 星 场 观测 得 到 的 扭曲 模型 下 接 用 于 稀 玖 星 图 像 的 扭曲 改正 。 

具体 的 拼接 步骤 如 图 1。 输 入 图 像 可 以 是 原始 的 FITS 图 像 ， 也 可 以 是 去 扭曲 后 ( 即 扭曲 改正 之 后 ) 
的 FITS 图 像 。 


输入 星 像 


1 1 底片 常数 模型 的 求解 TERE et |) eae 
对 于 每 帧 CCD RR, AR FAN H ZUR Ae i 
建立 标准 坐标 与 像素 坐标 的 关系 。 该 模型 中 包含 MAS 预 处 理 一 > ER 一 >| 灰 度 分 配 
了 平移 、 旋 转 、 放 缩 和 切 变 等 仿 射 变 换 的 影响 ， 
对 应 于 指向 的 不 准确 ， 微 小 大 气 折 射 较 差 改 正 ， 图 1 图 像 拼 接 软件 的 处 理 经 过 和 算法 
光 行 差 较 差 改正 等 影响 。 Fig. 1 Processing and algorithmic flow in image combination 


具体 地 ， 星 像 天 球 坐 标 (a, 5) 转换 为 标准 坐标 (E, n) 的 关系 方程 为 


cosdsin(a — A) 
sinDsinéd + cosDcosdcos(a — A) 


(1) 


_ cosDsinéd — sinDcosdcos(a — A) i 
sinDsinô + cosDcosdcos(a — A) 


FOP, (A, D) 是 投影 平面 与 天 球 切 点 的 天 球 坐 标 。 选 取 图 像 中 靠近 视 场 中 央 的 参考 星 理 论 位 置 作为 
切 点 坐标 。 
六 篆 数 模型 是 指标 准 坐标 (上 , n) 和 像素 坐标 (x, y) 的 关系 ， 关 系 方程 为 
E€=ax+by+c 
n=dx+ey +f” 
其 中 ，w、D、c、d、e、j 是 6 个 底片 常数 。 实 际 上 ， 采 用 最 小 二 乘 拟 合 求解 (2) 式 中 的 6 个 底片 常数 。 
1.2 像 系 坐标 转换 
选取 观测 天 区 中 靠近 中 心 部 分 的 图 像 作为 基准 图 像 ， 将 所 有 其 他 图 像 转换 到 基准 图 像 。 图 2 是 转 
换 过 程 的 流程 图 。 具 体 地 ， 对 待 转换 图 像 的 每 一 像素 点 的 位 置 ， 由 (2) 式 计算 其 标准 坐标 ， 由 (1) 式 
的 逆 运 算 再 转换 到 天 球 坐标 ， 进 而 从 天 球 坐标 转换 到 基准 图 像 中 ， 得 到 它们 在 基准 平面 的 像素 位 置 。 
其 中 ，(4, DD) 和 (4。, Do) 分 别 是 输入 图 像 的 切 点 和 输出 图 像 的 切 点 。 
1. ZK RE 和 个 
r 一 个 输入 像素 点 的 位 。 WAE > ea | 天球 村 
a 不 age ! ZAHID 式 (2) 转 换 “| 标准 坐标 | tia 
置 被 转换 成 基准 图 像 的 位 置 时 ， 它 的 像素 坐标 不 
一 定 是 整数 。 该 像素 的 灰 度 值 被 分 配 到 4 个 相 邻 
的 像素 中 。 如 果 感 光 融 的 填充 因子 为 100%， 可 | 输出 坐标 
以 按照 该 像素 在 基准 图 像 中 履 盖 不 同 像素 的 面积 


(2) 


输出 图 像 的 | 输出 切 点 为 (40, Do) 
ER 
ROAR PEK] 式 (D 转 换 


大 小 进行 灰 度 值 的 比例 分 配 。 通 过 上 述 计 算 ， 把 | A2 坐标 转换 过 程 
整 幅 输入 图 像 的 各 像素 的 灰 度 值 分 配给 输出 像 Fig.2 Procedure of coordination transformation 


系 ， 输 出 图 像 ( 即 基 准 图 像 ) 就 获得 了 该 图 像 的 灰 度 信息 。 继 续 处 理 下 一 幅 图 像 时 ， 在 输出 图 像 原 有 
的 灰 度 值 上 继续 累加 。 随 春 累 加 的 层 数 不 断 增加 ， 灰 度 值 也 越 来 越 大 。 为 了 归 一 化 ， 将 灰 度 值 除 以 拼 
接 图 像 的 帧 数 使 整个 输出 图 像 各 个 位 置 灰 度 值 转换 到 相同 的 曝光 水 平 。 
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2 观测 资料 


使 用 的 观测 资料 是 丽江 2.4m 望远镜 (观测 终端 为 YFOSC) 于 2013 年 2 月 5 日 拍摄 的 共计 58 IE 
动 观测 的 CCD 图 像 。 云 南 暗 天 体 相 机 及 分 光 仪 (Yunnan Faint Object Spectrograph and Camera, YFOSC) 
的 主要 参数 见 表 1。YFOSC 采用 背 照 式 CCD ， 该 CCD 结合 了 一 个 极 低 噪声 的 放大 器 。 
表 1 YFOSC 的 主要 技术 参数 
Table 1 Main technical parameters of YFOSC 


光学 系统 CCD RMI at 
准 直 镜 焦 距 . 222. 44 mm 有 效 靶 面 大 小 : 2 048 x 2 048 pixel? 
成 像 镜 焦距 : 114 mm 像 元 尺寸 ; 13.5x 13.5 pm? 
成 像 镜 焦 比 : /4. 1 每 个 像素 对 应 的 天 空 角 : 0. 285 arcsec 
焦 面 尺度 : 20. 96 arcsec/ mm 填充 因子 . 100% 


每 帧 图 像 的 曝光 时 间 均 为 30 s， 并 采用 B 滤 光 片 。 由 于 抖动 观测 ， 有 些 星 像 多 次 重复 出 现在 不 同 
图 像 中 ， 同 一 星 像 最 多 被 观测 51 次 。 做 平 场 处 理 并 裁剪 后 的 图 像 大 小 为 1 900 pixel x 1 900 pixel, %5 
9' x 9'。 用 文 [12] 的 方法 迭代 计算 得 到 受 扭曲 的 图 像 的 扭曲 量 进而 改正 得 到 不 受 扭曲 的 图 像 。 


3 ”结果 及 分 析 


图 3 左边 子 图 是 由 所 有 58 幅 不 受 扭 曲 的 输入 图 像 拼接 得 到 的 输出 图 像 ， 它 的 大 小 是 4 478 pixel x 
4 137 pixel， 周 围 空白 部 分 是 没有 拍摄 到 的 部 分 。 图 3 右边 子 图 是 左边 中 间 黑 色 和 矩形 部 分 的 放大 。 由 
于 Apophis 的 运动 轨迹 接近 于 一 条 直线 ， 并 且 拼 接 图 像 中 刚好 有 拍 成 一 列 的 多 个 星 像 ， 经 与 历 表 比 
BE, 证实 确 为 Apophis 的 轨迹 。 图 像 拼 接 相当 于 把 多 帧 图 像 的 信息 累加 来 增长 曝光 时 间 ， 这 种 增长 与 
重 琶 层 数 成 正比 。 考 虑 到 油 松 统计 规律 ， 图 像 的 拼接 相当 于 增加 了 曝光 时 间 ， 因 此 ， 信 咯 比 显著 提 
高 。 这 种 增加 信 噪 比 的 方法 避免 了 长 曝光 观测 导致 亮 星 饱和 溢出 。 如 图 4， 随机 选取 输出 的 拼接 图 像 
(H8) 的 一 块 区 域 与 原始 的 单 幅 输入 图 像 (左边 ) 的 同 区 域 作 比较 。 用 短 曝 光 担 摄 时 ， 星 像 非 常 模 糊 ， 有 时 
根本 无 法 识别 。 例 如 ， 输 入 图 像 ( 图 4 左边 ) 加 圆圈 部 分 看 似 只 存在 一 些 噪 声 , 但 输出 图 像 ( 图 4 右边) 


ER 4 
© 


图 3 不 受 扭曲 输出 图 像 及 其 局 部 放大 后 的 细节 
Fig. 3 The output image and the partial enlarged details 
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的 星 像 明显 清晰 可 见 。 通 过 测 光 发 现 ， 在 拼接 图 像 中 最 暗 的 星 像 星 等 估计 约 22 等 (B LSE). H 
SAOImage DSO 天 文 图 像 软件 绘制 该 星 像 在 不 同 输入 图 数量 的 拼接 图 像 中 的 光度 分 布 曲线 。 如 图 5， 


10 100 原始 单 幅 图 像 5 幅 拼接 成 像 


10 000 oon 


9 900 9 960 


9 800 


10 000 10 000 


9 980 


9 980 


9 960 9 960 


9 940 9 940 


图 5 不 同 输 入 帧 数 拼接 的 茶 暗 星 的 光度 分 布 曲 线 


Fig.5 Curve of intensity distribution of a faint under different input-frame numbers of combination 
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ER T 4 NRA X AT Te GRE ay A HAR, DUA ea TT BUEL, MERR X HAA hn ( AYE EB 17 
个 像素 ) 4 子 图 分 别 为 单 幅 原始 图 像 、5 Ba, 10 WE RRS APTA (58 幅 ) 图 像 拼 接 后 
的 灰 度 分 布 。 明 显 地 ， 重 登 层 数 越 多 ， 即 便 是 暗 星 的 形状 也 盘 近 高 斯 分 布 。 工 轴 方 向 的 光度 分 布 具有 
类 似 的 规律 。 

下 面 研究 图 像 扭曲 改正 对 拼接 图 像 的 影响 。 根 据 昨 表 UCAC4 对 每 一 幅 CCD 图 像 中 的 星 像 进行 归 
算 ， 可 将 输入 图 像 的 每 个 星 像 的 量度 坐标 换算 成 在 输出 图 像 中 的 量度 坐标 。 由 于 每 个 星 像 会 出 现在 不 
同 的 输入 图 像 中 ， 就 有 多 个 量度 坐标 ， 将 这 些 量 度 坐 标 换算 到 基准 平面 并 取 平 均值 作为 它 的 理论 位 
置 。 然 后 求 出 每 个 星 像 在 拼接 图 像 中 搜索 的 位 置 ， 并 计算 它 与 理论 位 置 之 差 。 在 x 方向 ( 赤 经 方向 ) 
My 方向 ( 赤 纬 方向 ) 的 搜索 位 置 与 理论 位 置 的 差分 别 简 记 为 Ax 和 Ay， 单 位 为 pixel, WEEN 
数 的 视 场 中 央 区 域 计 算 Ax 和 Ay， 实 验 结果 如 图 6， 除 了 个 别 暗 星 偏 得 比较 大 以 外 ， 其 余 偏差 都 很 
小 ，Ax 和 Ay 的 期 望 值 在 + 0. 02 pixel 之 内 。 对 于 受 扭 曲 的 拼接 (输入 图 像 是 受 扭曲 的 )，Ax 和 Ay 的 
标准 差分 别 是 0. 027 pixel( 约 8 mas) 和 0. 029 pixel。 而 对 于 去 扭曲 的 拼接 (输入 图 像 做 了 扭曲 改正 )， 
Ax 和 Ay 的 标准 差分 别 减 小 到 了 0. 014 pixel 和 0. 018 pixel。 可 见 扭曲 改正 的 效果 明显 。 当 然 ， 这 里 的 
误差 仅 代表 拼接 过 程 引 入 的 误差 ， 并 非 最 终 的 测量 精度 。 图 7 是 输入 图 像 与 拼接 图 像 的 星 像 位 置 残 差 
的 标准 差 (SD) 。 从 中 也 反映 出 扭曲 改正 对 拼接 结果 的 影响 显著 。 有 几 个 点 的 标准 差 侦 高 ， 这 是 输入 
图 像 间 的 某 些 重 僵 部 分 背景 差异 较 大 或 有 运动 日 标 穿 过 星 像 造 成 的 。 


* 
% Xo 


Ax/pixel 
Ay/pixel 


Age 


图 6 拼接 图 中 的 像 位 置 与 计算 位 置 的 偏差 


Fig.6 The deviation between the imaging position of combination and the calculated position 


在 和 X 轴 方向 的 SD/pixel 
在 立轴 方向 的 SD/pixel 


图 7 拼接 图 像 中 星 像 的 标准 差 


Fig.7 Standard deviation of stars in the combination 


2 期 阳 强 等 ; CCD 图 像 拼 接 试验 255 


4 4 论 


本 文 尝试 了 高 精度 CCD 图 像 的 拼接 。 试 验 发 现 ， 经 扭曲 改正 后 的 图 像 能 得 到 良好 质量 的 拼接 图 
像 。 具 体 表现 在 暗 星 数量 的 增多 和 信 噪 比 的 明显 改善 ,运动 目标 更 易 识别 。 但 计算 效率 需要 优化 改 
进 。 完 成 58 幅 图 的 拼接 时 ， 双 核 4 G 内 存 的 计算 机 运行 耗 时 长 达 10 多 分 钟 。 速 度 慢 是 因为 迭代 计算 
是 逐 像素 的 插值 整合 ， 每 次 整合 都 涉及 坐标 的 计算 和 插值 相交 面积 的 计算 ， 而 每 幅 图 又 有 数 百 万 的 像 
素 。 另 外 ， 将 继续 进行 背景 匹配 的 研究 ， 进 一 步 提高 精度 。 
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An Experiment of CCD Frame-Combination 
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Abstract: In order to obtain astronomic images with a wide field of view and high precision, a method of 
combining a few CCD (Charge Coupled Device) frames taken by a telescope is performed in this thesis. Peng 
et al.'s geometric distortion model is used to correct original images during the pre-processing phase. 
Coordinate transformation from original images to the final combined frame uses a six-parameter model together 
with the Catalog UCAC4. Unlike hardware image mosaic, the gray value of each pixel in the original frame 
needs to be redistributed to the pixels of a unique reference frame. The algorithms of this frame-combination 
are described in detail in this article. Furthermore, this method has been tested with many CCD frames taken 
by the 2. 4m Telescope at Yunnan Observatories. Results show that this method can yield a quite good CCD 
frame with a wide field of view. Specifically, moving objects can be detected easily; in terms of faint ones, 
there is significant improvement in their SNR (Signal to Noise Ratio) as well. Test results also show that the 
star positional errors introduced by the frame-combination are reduced to about 0. 02 pixels when raw CCD 
frames are corrected for their geometry distortions. This experiment therefore confirms the validity of Peng et 
al.'s geometric distortion solution and the necessity of correcting geometric distortions. 


Key words: Combination; Wide field of view; Geometry distortion; CCD frames 


